Effect of nonequilibrium phonons on hot-electron spin relaxation in n-type GaAs 

quantum wells 



P. Zhang and M. W. Wifl 



O 

(N 
C 

a 
in 



CO 



c 

o 
o 



> 

(N 

\o 

o 
o 



X 



Hefei National Laboratory for Physical Sciences at Microscale and Department of Physics, 
University of Science and Technology of China, Hefei, Anhui, 230026, China 
(Dated: January 15, 2010) 

We have studied the effect of nonequilibrium longitudinal optical phonons on hot-electron spin 
relaxation in ra-type GaAs quantum wells. The longitudinal optical phonons, due to the finite 
relaxation rate, are driven to nonequilibrium states by electrons under an in-plane electric field. The 
nonequilibrium phonons then in turn influence the electron spin relaxation properties via modifying 
the electron heating and drifting. The spin relaxation time is elongated due to the enhanced electron 
heating and thus the electron-phonon scattering in the presence of nonequilibrium phonons. The 
frequency of spin precession, which is roughly proportional to the electron drift velocity, can be 
either increased (at low electric field and/or high lattice temperature) or decreased (at high electric 
field and/or low lattice temperature). The nonequilibrium phonon effect is more pronounced when 
the electron density is high and the impurity density is low. 

PACS numbers: 72.25.Rb, 71.10.-w, 63.20.kd 



I. INTRODUCTION 

Understanding spin relaxation is an important issue for 
the possible application of spintronic devices.— ~— Among 
different kinds of spin relaxation mechanisms*^— scatter- 
ing plays an essential role. In general cases, phonons 
are assumed to form an equilibrium bath when carrier- 
phonon scattering is considered. This treatment works 
well when the carrier system is near the equilibrium. 
If the carriers are far away from the equilibrium (e.g., 
driven by an electric field or excited by a laser beam), 
phonons can be driven to run away from their equilib- 
rium states significantly by carriers when the carrier en- 
ergy relaxation mainly goes through the phonon emis- 
sions and the phonon relaxation time is comparable with 
(or longer than) the carrier-phonon scattering time. The 
nonequilibrium phonons in turn are able to affect the 
electron dynamics, including the spin relaxation. In fact, 
the hot-electron transport with nonequilibrium phonons 
has been investigated^— showing that the calculated 
electron energy loss rate and mobility fit better with ex- 
perimental data than those obtained with the equilibrium 
phononSi 10 : 13 : 14 These studies also indicate that it is nec- 
essary to treat phonons as nonequilibrium ones in the 
hot-carrier system, and the nonequilibrium phonons may 
affect spin relaxation via modifying the carrier heating 
and drifting. 

The hot-electron spin relaxation/dephasing has been 
studied theoretically in both (001) quantum-well 
structures^— and bulk materials^ by means of the ki- 
netic spin Bloch equation (KSBE) approach. 4 The spin 
relaxation/dephasing time is found to increase with elec- 
tric field when both the temperature and electric field are 
low, especially in high mobility samples*!^— When the 
electric field is high [for which the multi-subband (in con- 
fined nanostructures)iS and/or multi- valley^ effect have 
to be taken into account], the spin relaxation/dephasing 



time decreases with electric field^£ ~ 19 ' 21 In these studies 
the phonons are treated as equilibrium ones. This work is 
to investigate the influence of nonequilibrium phonons on 
hot-electron spin relaxation in an n-type GaAs quantum 
well, where the spin-orbit coupling term is the Dressel- 
haus typ o 19 ' 22 and the spin relaxation is limited by the 
D'yakonov-Perel' mechanism. 5 

The paper is organized as follows. In Sec. II we 
set up the model and the KSBEs with nonequilibrium 
phonons. In Sec. Ill the effect of nonequilibrium phonons 
on spin relaxation is investigated. Finally, we conclude 
in Sec. IV. 



II. MODEL AND KSBES 

We start our investigation from an n-type (001) ||z 
GaAs quantum well with an in-plane electric field. The 
well width a = 5 nm. Only the lowest subband is relevant 
with the proper electron density N e , lattice temperature 
Tl and electric field E. Due to the electron localiza- 
tion in the z-direction, the electron-phonon coupling is 
spatially inhomogeneous, i.e., the emission and absorp- 
tion of phonons mainly occur in the well where electrons 
have substantial density. If the phonon relaxation is fast 
enough or the phonons [particularly, the acoustic (AC) 
phonons] can easily penetrate through the well inter- 
faces, these phonons can be deemed as in equilibrium 
with the bulk modes. In our study we assume that the 
AC phonons keep in equilibrium and the longitudinal- 
optical (LO) phonons are nonequilibrium^ r 13 ' 1 5 In or- 
der to investigate the spin relaxation of electrons which 
are inhomogeneously coupled with the nonequilibrium 
LO phonons, we combine the rate equation of the LO 
phonons [Eq. (J2J], described as "quasi-2D" ^MLM. with 
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the electron KSBEs [Eq. ([T])]: 4 

dpk 







dt 


dt 




<9n q 


dt 


dt 



coh 



dpi, 
dt 



(1) 

(2) 



In Eq. (p}, Pk represent the density matrices of elec- 
trons with in-plane momentum k, whose diagonal terms 
Pkercr = /k<r (c = ±1/2) represent the electron distribu- 
tion functions and the off-diagonal ones /0 k i_i = pt_n 

2 2 k 2 2 

describe the inter-spin-band correlations for the spin co- 
eE • VkPk are the driving terms from 



herence. ^ 



dri 



are the coherent terms 
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the external electric field. 

describing the coherent spin precessions due to the effec- 
tive magnetic fields from the Dresselhaus ter m 19 i 22 and 
the Hartree-Fock Coulomb interaction, as well as the op- 
tional external magnetic field in the Voigt configuration. 

^r- stand for the scattering terms of electrons, in- 

scat 

eluding the electron-LO/AC phonon, electron- impurity 
and electron-electron Coulomb scatterings. n q in Eq. @ 
are the distributions of quasi-2D LO phonons with in- 



plane momentum q. 
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stand for the scattering 



terms of the LO phonons, including the phonon-phonon 
and phonon-electron scatterings. Expressions of the co- 
herent and scattering terms of electrons are given in 
detail in Refs. [I?] and [2(], except that the electron-LO 
phonon scattering term should be slightly modified here 
as the LO phonons are described as quasi-2D (this modi- 
fication makes no difference when the LO phonons are in 
equilibrium) i 10 i 13 ' 14 The electron-LO phonon scattering 
term in Eq. (fT]) reads 
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and <Sk(<; >) obtained by interchanging > and < from 
5k(>j<)- The scattering term in Eq. (J5|) reads 
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In the above equations p k 

riq = tiq and n q = n q + 1. £k — ti 2 ]< 2 /2m* rep- 
resents the energy of electron with momentum k and 
effective mass m* — 0.067mo, and hajQ = 35.4 meV 
is the energy of the LO phonons^ A is the area of 



the well layer. Af 2 = \ J2 Qz 9%, l^fe)| 2 is the effec- 
tive electron-LO phonon scattering matrix element with 

9qq. = 2tiie«$+qi) ^ ~ ^Q 1 )^ L iS th(3 SiZ(3 ° f the Sam " 

pie along the z-direction. kq — 12.9 and Koo — 10.8 are 
the relative static and high-frequency dielectric constants 
respectively^ and eo is the vacuum dielectric constant. 

stands for the form 
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factor under the infinite-depth well approximation. The 
first term on the right hand side of Eq. ([5]) represents the 
contribution from the phonon-phonon scattering in relax- 
ation time approximation. n q = [exp(fiwo/fc_B7L) — 1] _1 
is the number of quasi-2D LO phonons in equilibrium 
with the AC phonons. The population relaxation time 
T pp is contributed by anharmonic lattice vibrations (es- 
pecially the third-order anharmonicity) , which in prin- 
ciple depends on the phonon momentum and lattice 
temperature. Moreover, distinctly heated nonequilib- 
rium LO phonons (depending on the heating and relax- 
ation of electrons) may have different relaxation times. 
In spite of these intricate factors involved in the LO 
phonon relaxation, we assume t pp to be a constant only 
depending on the lattice temperature, by adopting r pp 
fitted by Vallee and Bogani from the time-resolved co- 
herent anti-Stokes Raman scattering experiment^ It 
gives t pp (T l ) = T° p /{1 + [cxp(0.2hu /k B T L ) - l]" 1 + 
[exp(0.8ftw /A: B T L ) - l]" 1 } with r p ° p = Tpp {0) « 9 ps. 24 
This formula in fact depicts the dominant decay route 
of an LO phonon near the center of the Brillouin zone 
into a transverse AC phonon and a different LO phonon 
at the L critical point of the Brillouin zone. The re- 
laxation time approximation with a constant t pp related 
only to the lattice temperature has been widely utilized 
in the study of hot-electron transport with the presence 
of nonequilibrium LO phonons 



III. RESULTS 

We numerically solve the KSBEs following the scheme 
mainly laid out in Ref. Il7l with the rate equation of 
the LO phonons discreted in the momentum space in 
a way similar to that for electrons. The impurity den- 
sity is set as zero and the electric field E = — £x with 
E > 0. No magnetic field is applied except other- 
wise specified. The initial conditions at time t = 
are chosen as the steady-state solution of Eqs. (flj and 
@ in the absence of the spin-orbit coupling in the co- 
herent term ^ 1 ^ Numerically, they are prepared 

coh 

from a state at t = —to (to > 0) with fka(-to) — 
{expp 2 k 2 /2m* - Po)jk B T L \ + l}" 1 , p krJ - a (-t ) = 



and 7i q (— to) = nl. Here p a are the electron chem- 



1 — p ki ical potentials determined by -j ^ k Tr[pk(— t )] = A e 



and i Sk Tr [Pk(-^o)^] 
the spin polarization and N, 



N e P Q , where Pq 



4 x 10 11 cm 



= 0.05 is 
" 2 is the 



electron area density. With the driving from the elec- 
tric field and the scattering, the system reaches a steady 
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state at time t = 0. After time t = 0, the spin-orbit 
coupling in the coherent term is switched on and elec- 
tron spins begin to relax with an initial spin polariza- 
tion P(t = 0) = Pq = 0.05. The spin relaxation time r 
is obtained from the time evolution of spin polarization 
-^W = ~Xn~ Tr[pk(i)c z ], the electron drift velocity 
is the steady value of v(t) = 4 ^ k Tr[pk(i)]^-k/m* = 
v x (t)St and the hot-electron temperature T e is fitted out 
from the Boltzmann tail of the steady-state electron 
distribution! 17 ' 19 



states by the hot-electrons. In the case with Tl — 200 K 
and E = 0.3 kV/cm [Fig. [Tib)], the LO phonons with 
q x > are emitted (and form a peak in the (^-positive 
momentum region) and those with q x < are absorbed 
(and form a valley in the g^-negative momentum region) . 
With the decrease of Tl and/or the increase of E, the val- 
ley in the ^-negative momentum region is suppressed or 
even disappears, as shown in Fig. QJd) f° r the case with 
Tl = 50 K and E = 1 kV/cm. In any case, the total 
phonon density increases and a net positive momentum 
is gained by the LO phonons. 
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FIG. 1. (Color online) Typical steady-state distributions of 
nonequilibrium electrons [(a) and (c)] and LO phonons [(b) 
and (d)] in momentum space, (a) and (b): T = 200 K and 
E = 0.3 kV/cm; (c) and (d): T = 50 K and E = 1 kV/cm. 
Note that in (b) and (d) the LO phonon distributions n q are 
rescaled by a factor 2 and =, respectively. 



In Fig. [TJ the typical steady-state distributions of the 
nonequilibrium electrons and LO phonons in momentum 
space are plotted. It is shown by Fig. [TJa) and (c) that 
under the electric field along the —x-direction, the elec- 
trons gain a drift velocity along the x-direction. Fig- 
ure QJb) and (d) show the corresponding distributions of 
nonequilibrium LO phonons. The LO phonons with ei- 
ther very large or small momenta [e.g., in the edge or 
center of the momentum space shown in Fig. QJb) and 
(d)] are in equilibrium with the AC phonons. It is seen 
from the figure that the equilibrium distributions of the 
LO phonons are different due to the distinct lattice tem- 
peratures [200 K for Fig. [TJb) and 50 K for Fig. [TJd)]. 
It is interesting to see that the LO phonons with medi- 
ate momenta are driven far away from their equilibrium 
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FIG. 2. (Color online) Electric field dependences of drift ve- 
locity v x (a), hot-electron temperature T e (c) and spin relax- 
ation time r (e), calculated with the equilibrium LO phonons. 
The ratios of the quantities v x , T e and r with the nonequilib- 
rium LO phonons to those with the equilibrium ones, r) Vx , rfr e 
and r\ T , are shown in (b), (d) and (f), respectively. In (e), the 
two curves with open squares and open circles are the elec- 
tric field dependences of r with a magnetic field (B — 2 T) 
along the —x-direction for Tl = 50 K and 300 K, respec- 
tively [at T L = 50 K (T L = 300 K), when E > 0.05 kV/cm 
(E > 0.3 kV/cm) the spin relaxation time coincides with that 
without magnetic field and thus not shown]. 
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In Fig. [2a) the electric field dependence of electron 
drift velocity v x at various lattice temperatures is plot- 
ted, with the LO phonons treated as the equilibrium ones 
in the calculation. The ratio of v x obtained with the 
nonequilibrium LO phonons to that with the equilibrium 
ones is shown in Fig. [2jb) . From Fig. [2jb), one notices 
that when the nonequilibrium phonon effect is taken into 
account, v x can be either increased or decreased. In fact, 
the influence of nonequilibrium phonons on hot-electron 
transport consists of two competing effects: the reab- 
sorption of momentum from the nonequilibrium phonons 
tends to increase the electron mobility, while the en- 
hanced electron heating strengthens the electron-phonon 
scattering (including both the electron- AC phonon and 
electron-LO phonon scatterings) and tends to decrease 
the electron mobility^ ' 11 ' 13 Generally the former (latter) 
effect dominates in the regime with low (high) electric 
filed and/or high (low) temperature^ ' 11 ' 13 as indicated 
in Fig. HJb) • As mentioned previously, when the electric 
field is low and/or the lattice temperature is high, the val- 
ley of the LO phonon distribution in ^-negative momen- 
tum space is pronounced, thus it substantially suppresses 
the back scattering of electrons in momentum space by 
absorbing the q x -negative phonons and hence makes the 
electron distribution more forward-peaked. Finally, when 
the lattice temperature is high enough (e.g., Tl = 300 K), 
the effect of nonequilibrium phonons becomes weak and 
less sensitive to the electric field, mainly due to the 
shorter LO phonon relaxation time (t pp rs f .9 ps when 
T L = 300 K)i2± 

From Fig. [He), where the electric field dependence 
of the hot-electron temperature T e under various lattice 
temperatures is shown, one indeed finds that the heating 
of electrons by the electric field is quite obvious when 
the electric field is high and the lattice temperature is 
low . 17 ' 19 ' 21 From Fig. [U(d), where the ratio of T e ob- 
tained with the nonequilibrium LO phonons to that with 
the equilibrium ones is plotted, one finds that with the 
nonequilibrium phonon effect considered, electrons are 
further heated as expected.— ~— Moreover, when Tl is 
low, r]T B shows a nonmonotonic behavior. That is caused 
by the decay of the heating efficiency with the increase of 
electric field in the presence of nonequilibrium phonons: 
With the increase of electric field, the number of the 
LO phonons increases and the rate of electron energy 
relaxation through the electron-LO phonon scattering in- 
creases as well. This effect is more pronounced when the 
lattice temperature is low, where the nonequilibrium LO 
phonons can be considerably accumulated with the in- 
crease of electric field, due to the long phonon relaxation 
time (t pp — 7.3 ps when Tl = 50 K)2£ as well as the 
small equilibrium phonon distribution. 

The electric field dependence of the spin relaxation 
time r under various temperatures is plotted in Fig.|2Je). 
From the figure, one notices that r generally increases 
with E in the regime under investigation. Two rea- 
sons lead to this phenomenon: (I) Under the electric 
field along the —x-direction, a net effective magnetic 



field along the —x-direction is induced via the Dressel- 
haus spin-orbit coupling . 17 ' 19 With this effective mag- 
netic field, spins begin to precess around it and thus the 
in-plane spin relaxation is mixed with the out-of-planc 
onei 25 ' 26 The two-dimensional electron system in (001) 
GaAs quantum well has an in-plane spin relaxation rate 
smaller than the out-of-plane one in the framework of 
the D'yakonov-Perel' relaxation mechanism. Thus when 
the electric field is applied, r is increased due to the ef- 
fective magnetic field . 25 ' 26 The effective magnetic field 
decreases with the increase of Tl, because it is propor- 
tional to v x 17 i 19 and v x decreases with increasing Tl [as 
shown in Fig. [2(a)]. Therefore the mixing of the in-plane 
and out-of-plane spin relaxations is obvious in the low 
temperature regime. In fact, when E = 0.05 kV/cm, the 
effective magnetic field is ~ 1 T when Tl ^50-100 K and 
- 0.1 T when T L = 300 K. As a result, the effective 
magnetic field causes an abrupt increase of r [Ref. [25| 
with the increase of E from to 0.05 kV/cm for the 
cases with Tl ^50-100 K but a slow increase of r with 
the increase of E from to 0.3 kV/cm for the case with 
Tl = 300 K. (II) The heating of electrons by the electric 
field enhances the electron-phonon scattering and thus 
increases the spin relaxation time r in the strong scat- 
tering limit.— li 7 - - — This effect, only important in the low 
temperature regime where the heating effect is strong 
[as shown in Fig. [DJc)], is responsible for the continuing 
increase of t with E when E > 0.05 kV/cm and Tl ^50- 
100 K. To make the underlying physics depicted above 
more pronounced, a magnetic field B — 2 T is applied 
along the —x-direction for the cases with Tl — 50 K and 
300 K. The corresponding electric field dependences of r 
are plotted as curves with open squares (Tl — 50 K) and 
open circles (T L = 300 K) in Fig.^e). With this large ex- 
ternal magnetic field, the in-plane and out-of-plane spin 
relaxations are efficiently mixed even when E = 0. One 
finds that for the case with Tl = 300 K, r almost keeps 
unchanged with the increase of electric field, while for the 
case with Tl = 50 K r keeps on increasing with E due to 
the strong heating effect. Finally, it is noted that when 
Tl = 300 K, there is a marginally decreasing tendency 
of t with E when E is near 1 kV/cm. This is caused by 
the enhanced inhomogencous broadening of the effective 
magnetic field from the Dresselhaus spin-orbit coupling 
due to the drifting and heating of the electric field^ 7 - - — 
This effect is easier to take place when the lattice tem- 
perature is high i 17 ' 19 

In Fig. [2{f), the ratio of the spin relaxation time ob- 
tained with the nonequilibrium LO phonons to that with 
the equilibrium ones is shown. With the nonequilib- 
rium LO phonons taken into account, t is generally in- 
creased due to the strengthened electron-phonon scat- 
tering. Therefore the increase of r corresponds to the 
increase of T e , as shown in Fig.[^d). However, when the 
lattice temperature is high enough (e.g., Tl — 300 K), 
the modification on r induced by the nonequilibrium 
phonons can not be seen. Moreover, as the spin preces- 
sion frequency is proportional to v x ( Refs . [l7l and [lit) and 
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v x is affected largely by the nonequilibrium LO phonons, 
the spin precession frequency has a modification with the 
magnitude roughly proportional to the change of v x . In 
Fig-EJa) we show the typical spin precession signals with 
the equilibrium and nonequilibrium LO phonons, respec- 
tively. In Fig. [3jb) , the typical ratio of spin precession 
frequency obtained with the nonequilibrium LO phonons 
to that with the equilibrium ones is plotted in the region 
of the electric field where v x is large enough that spin 
precession signal with clear periods can be distinguished. 
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FIG. 3. (Color online) (a) Fypical spin precession signals cal- 
culated with the equilibrium and nonequilibrium LO phonons. 
Note that the absolute value of spin polarization \P\ is in the 
logscale. T L = 100 K and E = 0.3 kV/cm. (b) The ratio 
of spin precession frequency obtained with the nonequilib- 
rium LO phonons to that with the equilibrium ones, shown 
at three typical temperatures in the electric field region where 
the spin precession signal can be clearly distinguished. 

Finally, further calculations show that the effect of 
nonequilibrium phonons on electron mobility, electron 
heating and electron spin relaxation decays with the 
decrease of electron density and the increase of impu- 
rity density. This can be easily understood as the LO 
phonons are driven to the nonequilibrium states by elec- 
trons and the increase of the electron-impurity scattering 
suppresses the effect caused by the electron-LO phonon 
scattering. 



rt-type (001) GaAs quantum wells. Under an in-plane 
electric field, the LO phonons can be driven away from 
their equilibrium states by electrons and then in turn af- 
fect the electron transport, electron heating as well as 
electron spin relaxation. 

In the presence of the nonequilibrium LO phonons, 
the electron drift velocity under electric field can be ei- 
ther increased (at low electric field and/or high lattice 
temperature) or decreased (at high electric field and/or 
low lattice temperature). This phenomenon is caused 
by the two competing effects: the momentum reabsorp- 
tion from phonons and the strengthened electron-phonon 
scattering . V 1 ? 13 The former tends to increase the elec- 
tron mobility whereas the latter tends to suppress it. The 
nonequilibrium LO phonons also impede the energy re- 
laxation of electrons and thus the electrons are further 
heated, especially when the lattice temperature is low. 
The nonequilibrium LO phonons effectively affect the 
hot-electron spin relaxation through the strengthening of 
the electron-phonon scattering, which tends to increase 
the spin relaxation time. This effect also dominates in 
the low temperature regime. Moreover, as the spin pre- 
cession frequency under the electric field is proportional 
to the electron drift velocity, it can be either increased 
or decreased when the nonequilibrium LO phonons are 
taken into account. Finally, it should be noticed that 
the effect of nonequilibrium phonons is more pronounced 
in systems with high electron density and low impurity 
density. 
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IV. COUCLUSION 



In this work, we have studied the effect of nonequi- 
librium LO phonons on hot-electron spin relaxation in 
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